Scaling of quantum Zeno dynamics in thermodynamic systems 
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We study the quantum Zeno effect (QZE) in two many-body systems, namely the one-dimensional 
transverse-field Ising model and the Lipkin-Meshkov-Glick (LMG) model, coupled to a central qubit. 
Our result shows that in order to observe QZE in the Ising model, the frequency of the projective 
measurement should be of comparable order to that of the system sizes. The same criterion also 
holds in the symmetry broken phase of the LMG model while in the model's polarized phase, the 
QZE can be easily observed. 
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I. INTRODUCTION 

Quantum Zeno effect (QZE) refers to the inhabitation 
of the unitary time evolution of a quantum system by re- 
peated frequent measurements [lH7[. It is a phenomena 
which is intrinsically related to the projection postulate 
in quantum mechanics. For an isolated quantum system, 
its state vector undergoes an unitary evolution according 
to the Schrodinger's equation. At the very beginning, 
the survival probability for the system to remain in the 
initial state changes quadratically with the elapsed time, 
while it decreases exponentially in a large time scale. If 
a measurement is performed to check whether the sys- 
tem is still in the initial state or not after some time, the 
projection postulate states that the system's state will 
collapse to the initial state immediately after the mea- 
surement. Therefore, when successive measurements are 
performed at a time interval within the quadratic decay- 
ing region, the unitary evolution of the system is found to 
be suppressed. In the limit of continuous measurement, 
the dynamic evolution of the system is frozen. 

The possibility of such an effect to be found in an un- 
stable quantum system was first pointed out by Fonda 



et al. pjj and Degasperis et al. [12J. They suggested 
that the decay rate of an unstable system may be de- 
pendent on the frequency of intermediate measurements. 
Later, Misra et al., who gave the effect its name, pre- 
sented a general formulism concerning the semigroups to 
analyze the QZE [l|. Experimentally, the QZE was first 
manifested in an radio-frequency transition between two 
laser-cooled Beryllium ion ground-state levels |13J |. 

In the recent decade, QZE has been studied intensively 
within the content of quantum optics. Among those anal- 
yses, the systems under consideration are only of a few 
levels in which the QZE can be easily observed by fre- 
quent measurements. However, little attention of the 
possibility to observe QZE in quantum many-body sys- 
tems has been paid so far. The primary motivation of 



our work is to investigate the QZE from the viewpoint 
of condensed matter physics and obtain criteria on how 
frequent the measurements should be compared to the 
system size in order to preserve the QZE in many-body 
systems. 



Mathematically speaking, we are interested in the scal- 
ing behavior of the short time limit of the survival prob- 
ability of the system's initial state under the influence 
of an external perturbation. We presented a general 
analysis and find that the leading term of the survival 
probability is equal to the fluctuation in the interaction 
Hamiltonian. Specifically, we take the one-dimensional 
transverse-field Ising model and the Lipkin-Meshkov- 
Glick (LMG) model as examples. The leading term of 
the survival probability in these two models are calcu- 
lated analytically and their scaling behaviors is obtained. 
In our analysis, we allowed these many-body models to 
couple to a central qubit. The significance of the cen- 
tral qubit is just to provide an external perturbation in 
the Hamiltonian and allow us to study how the system's 
evolution is affected by the external perturbation in the 
short time regime. A point to note here is that the central 
qubit only serve as a source of external perturbation and 
how it evolves is not of our interest. Moreover, we would 
like to point out that the leading term of the survival 
probability in our study is just the same as the linear 
response in the study of quench dynamics |8l4l0|| . 



'Electronic address: sjgu@phy.cuhk.edu. hk 



The paper is organized as follows: In Section UH we 
first present a general formulism for obtaining the sur- 
vival probability of the initial state in the short-time 
limit. Then in Section Mil and IIV1 we take the one- 
dimensional transverse-field Ising model and the Lipkin- 
Meshkov-Glick (LMG) model as examples respectively to 
illustrate the criteria for observing the QZE, in terms of 
the size dependence of the leading term of the survival 
probability in the short-time limit. Our analysis shows 
that in order to observe the QZE, the number of measure- 
ments have to be comparable to the size of the system. 
Finally, a summary would be given in Section [Vj 



II. FORMULISM 

Consider a quantum system interacting with the en- 
vironment, the Hamiltonian of the whole setup can be 
generally written as 



H = H a 



■SHj, 



(1) 



where Hq describes the initial Hamiltonian of the sys- 
tem and the environment such that Ho\m) = E m \m). 
Here {|m)} is a set of orthogonal basis constructed from 
the tensor product of the eigenstate of the system and 
that of the environment and E m is the corresponding 
eigenenergy. Hj is the coupling between the system and 
the environment with a small parameter 5 denoting the 
coupling strength. 

Suppose the setup is initially in a state \ip(0)} = |to}- 
It is allowed to evolve freely for a time r. According to 
the Schrodinger's equation, the state at time r is given 
by \iP(t)) = e~ lHT \m) and the probability of finding the 
setup in the initial state, i.e. the survival probability, is 
given by 



P(r) = b(T)| 2 H(<K0)kHr))| 2 



{m\e- l{Ho+SH,)T \m) 



(2) 



For small S, one can expand p(r) around 8 = and keep 
terms up to the second order. We have 
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where H 7 /" 1 = (n\Hi\m). Keeping terms to the lowest 
order of 5, Eq. ([2]) becomes 

P(t) ^l-25 2 J2 l ~ tJ^V E 2 l)T] \ H T n W (6) 



n^rn 



{E n - E m ) 2 



which is just the perturbative form of the Loschmidt echo 
[ljj. In the short time limit, Eq. © gives 



P(t) « 1 - S 2 t 2 J2 \ H ? 



(7) 



Now, suppose successive measurements are performed 
at every time interval r, and the system is allowed to 
evolve freely in-between consecutive measurements. The 
probability of finding the system in the initial state after 
a finite duration At = Qt, where Q is the number of 
measurements, is 



P(At) = \P(t 



\Q 



1 



*wi, 



where 



E 



H^ n \ 2 = {m\HJ\m) 



— (m\Hj\m) 



(8) 



(9) 



which is the fluctuation in the interaction Hamiltonian. 
Note that the projection postulate has been implanted in 
the first equality in Eq. ©. From Eq. ([5]), one realized 
that there exists a competition between \ arL d Q. For 
a few-level system, since x is finite, we have P(At) —t 1 
under the case of continuous measurements (Q — > oo). 
In other words, the system would remain in the initial 
state if the measurements are performed continuously. 
However, for a quantum many-body system, \ may not 
be finite but goes with some power of the system size. 
In the thermodynamic limit, \ is also infinite and the 
second term in Eq. (|8]) could not be simply ignored even 
if Q — > oo. The key motivation of our work is to obtain 
the scaling behavior of x in quantum many-body systems. 
From there, one can then predict the behavior of x i n the 
thermodynamic limit and draw criterions on how large Q 
should be in comparison to the system size in order to 
observe the QZE in the system. 

In the following, we take the one-dimensional 
transverse-field Ising model and the LMG model coupled 
to a central qubit as examples. Using Eq. ([S]), we calcu- 
lated x explicitly and extracted its size dependence. One 
can then determine the criteria on how large Q should 
be in order to observe the QZE in the models. 



III. ONE-DIMENSIONAL TRANSVERSE-FIELD 
ISING MODEL 

The Hamiltonian of the one-dimensional transverse- 
field Ising model reads 



H, 



Ising 



>?a?+i+ Aof), 



(10) 



where cr?, with q = x,y,z, is the Pauli matrix. It de- 
scribes a chain of spins interacting with the nearest neigh- 
bors in the z direction, while all the spins are subjected 
to a transverse external field with strength A in the x 
direction. 

The Hamiltonian in Eq. (ITU1) can be diagonalized ex- 
actly following three standard transformations [15|]: i.e. 
Jordan- Wigner transformation, Fourier transformation, 
and Bogoliubov transformation. Then the Hamiltonian in 
Eq. (|10p can be written into an diagonal form of -ffising = 
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FIG. 1: x as a function of A for the one-dimensional 
transverse-field Ising model coupled to a central qubit for 
system size N = 10, 16 and 20. The symbols denote data 
obtained from the numerical exact diagonalization while the 
lines are obtained from the analytical expression of \ in Eq. 
(|14fl . The analytical results agree with the numerical data. 



J2k £k (dkdk~- 1/2), where £fc = 1s/\ + A 2 — 2A cos k is the 
single quasi-particle energy. d k and d k is the fermionic 
annihilation and creation operator respectively. 

Now, consider a central qubit interacting transversely 
with the Ising model. The corresponding Hamiltonian 
reads 

H = Hiring -*^|e)<e|of, 

i 

= H Ising | g )(.g| + ( J ff Ising + ^ / )|e)(e|, (11) 



where -ffising is given by Eq. (fT0|) . and Hi = — Yli a i- 
\g) and |e) denotes the ground state and the excited state 
of the central qubit, respectively. The Loschmidt echo in 
this system has been investigated by Quan and his col- 
labrators [l6(. They found that the decay of Loschmidt 
echo is greatly enhanced at the critical point of the Ising 
model. In the following, we will follow similar approach 
as used by Quan et al., but focusing on the short-time 
limit of the Loschmidt echo. 

Without loss of generality, assume the Ising model 
is initially in some state |^(0)) and the central qubit 
is in a superposition state c g \g) + c e \e), where the co- 
efficients satisfie |c g | 2 + |c e | 2 = 1. The initial state 
of the whole setup is thus |-0(O)) = |0(O)) ® (c g \g) + 
c e \e}). The evolution of the Ising chain splits into two 
branches and the state vector of the setup at a later 
time r is given by \iJ>(t)) = c g e - * fflBta * T |0(O)) ® \g) + 
c e e- it - Hl ^+ 5H ^ T \(j)(0)) <g> |e). 

For simplicity, suppose the Ising chain is initially in the 
ground state \G), and the central qubit is in the excited 
state, i.e. c g = and c e = 1. The survival probability of 
the system as given in Eq. @ becomes 



P(r) = \(G\e 



-i(Hi s i ns +SH!)T 



\G)f 



(12) 



FIG. 2: The normalized \ as a function of A for the one- 
dimensional transverse-field Ising model coupled to a central 
qubit. The curve for different system sizes collapses to a single 
one. 



Note that in the eigen-basis of the diagonalized Hamilto- 
nian jrifc) such that d k dk\nk) = nk\nk), where nk = or 
1, the ground state of the Ising model can be written as 
\G)=I\ k> o\^0-k). 

By performing the Jordan- Wigner transformation, 
Fourier transformation, and the Bogoliubov transforma- 
tion, Hi can also be transformed into the eigen-basis 
{\n k )} as 

Hi = -N (13) 

+ y^ \2(cos e k )d\d k + i(sin6> fe )(dj c dt A , - d- k d k ) , 
fc 

where 9 k is the Bogoliubov angle satisfying tan(9fc = 
sinfc(cosfc — A) -1 . Using Eq. ©, we obtain 



fe>0 



4 sin 2 k 



A 2 — 2A cos k 



(14) 



Next, we would like to investigate the size dependence 
of x- Assume N is even, we introduce the density of 
states in A:-space in the limit of TV — > oo. For A = 1, one 
gets 



N f("-W N si n 2 fc 



X = — 

7T 



tt/JV 



1 — cos k 



■dk - N. 



(15) 



For A ^ 1, let < x < 1/2 and < A; = 2ttx < tt, then 
Eq. (fl4l becomes 



X = 2N 



sin 27ra; 



-dx. 



(16) 



ii 



1 + A 2 - 2Acos27ra;' 

Let z = e l2 ^ x , so sin(27r:r) = -^(z — -) and cos(27rx) = 
\{z+ i), the integral becomes the contour integral along 
a unit circle on the complex plane: 

N f (z 2 -l) 2 



X 



4ni 



z 2 (z-\)(\z-l) 



dz. 



(17) 



which can be easily evaluated using the residue theorem. 
One then obtains 



N/X 2 ~ N 

N ~ N 



for A > 1 
for A < 1. 



(18) 



To check the validity of the result obtained above, one 
may consider the limiting case where A — > 0. In this 
limit, the ground state of the system is a superposition 
of two polarized wavefunctions, i.e. -4=(| tt '"' t) + 
| J4- • • • 1)). Noting that the interaction Hamiltonain Hi 
only flips a particular spin at site i. This results in a 
wavefunction orthogonal to the original one. The second 
term in Eq. © is thus zero. For the first term, one can 
also easily verify that gives N using the property of the 
Pauli's matrices. Therefore, one obtains x = A^ as A ^ 0. 
On the other hand, form the analytical expression in Eq. 
(fT4"| . we also obtain the same result for A = as shown 
in Eq. (fl8j). 

Furthermore, we also perform the numerical exact di- 
agonalization and calculated \ using Eq. ([9]). Fig. Q] 
shows a plot of x as a function A for N — 10, 16 and 
20. The lines obtain from the analytical expression of \ 
agree with the numerical data (represented by symbols). 

Fig. [5] shows a plot of \ given by the analytical ex- 
pression in Eq. (|14p normalized over the system size as 
a function of A. The curve for different system sizes col- 
lapse to a single one, i.e. \ ~ N. This echoes the result 
obtained from the scaling analysis above. 

Returning to Eq. ©, we see that in order to observe 
QZE in the Ising model, Q 3> £N where £ is a constant 
independent of the system size. In the thermodynamic 
limit, in which N is of order 10 23 , Q has to be at least of 
this order or else the QZE breaks down. 

From Fig. [5J one can also realize that the normalized x 
changes manifestly across the point at A = 1, which is the 
critical point where a quantum phase transition(QPT) of 
the Ising model takes place [15[ . Differentiating Eq. (|T4|) 
with respect to A, one obtains 



dX 



E 

fc>0 



8(cosfc — A) sin 2 k 
(l + A 2 -2Acosfc) 2 ' 



(19) 



A plot of Eq. (|19|) normalized over the system size as a 
function of A in the vicinity of the critical point is shown 
in Fig. |31 d(x/N)/dX changes from positive to negative 
as the system goes across the critical point from below. 
This change becomes sharper and sharper as the system 
size increases. Fig. @] shows the second derivative of the 
normalized % as a function of A in the vicinity of the 
critical point. One can see that d 2 (x/N)/dX 2 exhibits a 
minimum at the critical point. The magnitude of this 
minimum is in fact scales with N, as shown in the inset 

of Fig. m 

The size dependence of d 2 {x/N)/dX 2 can also be ob- 
tained analytically. Differentiate Eq. (|T9"]) with respect 
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FIG. 3: The first derivative of the normalized x with respect 
to A as a function of A near the vicinity of the critical point for 
the one-dimensional transverse-field Ising model coupled to a 
central qubit. d(x/N)/d\ changes from positive to negative 
across the critical point of the Ising model. The change is 
faster for a larger system size. The curve for different system 
sizes coincide exactly at A = 1. 
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FIG. 4: The second derivative of the normalized \ with re- 
spect to A as a function of A near the vicinity of the critical 
point for the one-dimensional transverse-field Ising model cou- 
pled to a central qubit. The inset shows the scaling behavior 
of the minimum of d 2 (x/N)/dX 2 . The slope of the straight 
line is 1.007 ± 0.001. 



to A, one obtains 
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dX 2 
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FIG. 5: x as a function of A for the LMG model coupled to 
a central qubit with system sizes N = 800 and 1000 (7 = 
0). The dots denote data obtained from the numerical exact 
diagonalization while the smooth lines are obtained from the 
analytical expression of \ m Eq. (|38[) and (|40[) . The analytical 
results agree with the numerical data except for points close 
to A = 1. 



For A = 1, 

d 2 X 
d\ 2 



E 
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N 
N 2 



4 sin k 
1 — cos k 

(N-1)tt/N 
tt/N 



2 siri k 



(1 

2(1 



cosfc) 2 
cosfc) - 



1 + cos k 



1 — cos k 



dk, 

(21) 



We then have d 2 (x/N)/dX 2 ~ TV which diverges in the 
thermodynamic limit. 

Furthermore, from Fig. O we see that there exists a 
fixed point at A = 1 where d(x/N)/d\ for various system 
size coincide. This can also be shown analytically from 
Eq. dTSJ. For A = 1, Eq. (Q]5) gives 



d(x) 
dX 



Et 



2 sin 2 k 



cosk 

N r(N-l)n/N 

n Jn/N 
-N, 



(1 + cosk)dk 



(22) 



in which the average over system size gives —1. This 
fixed point reflects some kind of symmetry is hidden in 
the model and further exploration would be of interest. 



IV. LIPKIN-MESHKOV-GLICK MODEL 

The Hamiltonian of the LMG model reads 
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FIG. 6: (a) x as a function of A (b) The normalized x as a 
function of A for different system sizes for the LMG model 
coupled to a central qubit (7 = 0.5). The normalized x f° r 
various system sizes collapse to a single curve. 



The Hamiltonian describes a cluster of spins mutually 
interacting with each other in the xy plane. 7 is the 
parameter indicating the anisotropy in the interaction 
of the model. All the spins are embedded in an exter- 
nal field along the z direction. The prefactor 1/iV is to 
ensure a finite energy per spin in the thermodynamic 
limit. Introducing S K — J^ erf /2, where K = x,y, z, and 
S x = (S+ + S-)/2 and S y = {S+ - S-)/2i, Eq. ((53) can 
be rewritten as 



Hlmg 



-2XS, 
1 



5^(1 -7) (si 



s 2 J 



N 



N (i + i){s 2 -s 2 - 2 



(24) 



l<3 



In the large N limit and considering the low excitation 
spectrum, the above Hamiltonian can be diagonalized 
through a semi-classical approach as in the following jL7l| : 

(i) Perform a rotation of the spin operators around the 
y axis in order to bring the z axis along the semiclassical 



magnetization. This is done by 



cos a sin a 

1 
— sin a cos a 





(25) 



where a — for A > 1 so that S = S, and a = cos -1 A 
for A < 1. 

(ii) Holstein-Primakoff transformation which trans- 
forms the spin operators into bosonic operators by 

S z =S-aU = N/2 - at a, (26) 

S+ =(2S-a^a) 1 / 2 a = VN{l-a^a/N) 1 / 2 a, (27) 
S_ = a^(2S-a^a) 1 / 2 = VNa^(l-a^a/N) 1 / 2 , (28) 

where a and a^ are bosonic annihilation and creation op- 
erators satisfying the commutation relation [a,cr] = 1. 
The second equality in above equations hold as a result 
of [H, S 2 } = and only the low energy spectrum is con- 
sidered. 

(hi) Bog oliubov transformation: After the above trans- 
formations, we arrive at a quadratic Hamiltonian and one 
can easily diagonalize it using the Bogoliubov transfor- 
mation: 



a = sinh ( - ) b^ + cosh I - ) 6, 



a 1 * = cosh I - I &t + sinh - ] b, 



(29) 
(30) 



where b and 6^ are also bosonic operators satisfying the 
same commutation relation as that of a and a' . Following 
the above transformations, one could finally diagonalize 
the Hamiltonian in Eq. (j2~4")l into -Hlmg = e oN + (\ — 
A/2 + VA 2 -4r 2 (&t& + 1/2), where for A < 1, 

1 + A 2 



eo = 



ei 



i-S 



A = 2 -A 2 -7 
. _ 7-A 2 



(31) 



For A > 1, 



= -A 

= 

= 2A - 1 - 7 ■ 

r =-i^r 

2 
The diagonalization condition is given by 

2r 

tanhw = — — , 

A 

from which we also have 

A 



(32) 



cosh( 



VA 2 - 4r 2 ■ 



sinhfl 



-2T 



VA 2 - 4r 2 



(33) 



(34) 



Similar to the case of the Ising model in the previous 
section, let's consider a central qubit interacting trans- 
versely with the LMG model. The corresponding Hamil- 
tonian reads 



H = H 



LMG 



*E 



e)(e\at 



(35) 



where H^mg is given by Eq. (|2"3"|) . Again, suppose the 
LMG model is initially in the ground state \G), and the 
central qubit is in the excited state. Following the same 
argument as in the Ising model, the survival probabil- 
ity of the system remain in the initial state in Eq. ^ 
becomes 



P(t) = |(G|e-^ LMG+ "^ T |G)| 2 , 



(36) 



where Hi = — J2i °~i- Similarly, in the eigen-basis {\n}} 
such that b*b\n) = n\n), where now n = 0, 1,2- ••, 
note that ground state of the LMG model is just \G) = 
|0). Following the diagonalization procedure mentioned 
above, Hi can also be expressed in the eigen-basis {|n)}. 
For A > 1, 

Hi = N + l-cosh8-2cosh(6)tfb 

-sinh(0)(6 t2 +6 2 ). (37) 

Using Eq. ©, (|33]). and (|37]). one obtains 

X = 2sinh 2 6>- N°. (38) 

For A < 1, 

Hi = -X(N - cosh(6>) + 1) 

+2 cosh(6»)6 t 6 + A sinh(6») (fe t2 + b 2 ) (39) 

+ y/N(l - A 2 )(& + + 6)(cosh(6»/2) + sinh(f9/2)). 
Using Eq. ©, ([HI), and (JUJl, one has 
X = 2A sinh 2 9 



+A^A(1 - A 2 )sinh6»(Vcoshf9 + 1 + \/cosh(9 - 1) 
+ (1 - A 2 ) N (cosh 9 + sinh 6>) 
- N. (40) 

Again, to check the validity of the analytical result, 
let's first consider the limiting case where 7 = and A — > 
0. In this limit, the Hamiltonian as given by Eq. (|2"3]) is 
similar to that of the Ising model except the interaction 
is over all spins here. Following similar argument as that 
for the Ising model, \ obtained from Eq. (|9]) also gives 
N. On the other hand, from Eq. (13"TT) and (|3"4")l . one can 
easily see that Eq. (J40I) also gives N. The two results 
are consistent. 

Moreover, the analytical result is also compared with 
the data from the numerical exact diagonalization. As 
shown in Fig. [5] the two data agree with each other 
except for A close to 1. This invalidity of the analytical 
result for finite system sizes near the critical point is due 
to the fact that we have only kept terms up to the lowest 
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FIG. 7: First derivative of the normalized x with respect to 
A as a function of A in the vicinity of the critical point for the 
LMG model coupled to a central qubit (7 = 0.5). The inset 
shows the scaling behavior of the minimum of d(x/N)/d\. 
The slope of the straight line gives 0.324 ± 0.001. 



FIG. 9: First derivative of the normalized \ with respect to 
A as a function of A in the vicinity of the critical point for the 
LMG model coupled to a central qubit (7 = 0). The inset 
shows the scaling behavior of the minimum of d(x/N)/d\. 
The slope of the straight line gives 0.328 ± 0.001. 
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FIG. 8: (a) x 3 ) The normalized \ f° r different system sizes 
as a function of A for the LMG model coupled to a central 
qubit (7 = 0). The normalized \ for various system sizes 
collapse to a single curve. 



order in the 1/N expansion of the HP transformation. To 
obtain more accurate results near the critical point, one 
have to go to higher order terms in the 1/N expansion 
of HP transformation [F7| or through the analysis of the 
Majorana polynomial roots [18| , In this report, our focus 
is on the scaling behavior of \ within the same phase 
of the model. The analytical analysis presented so far 
remains valid as long as one stays in the same phase. 

Fig. [6] shows a plot of x obtained from exact diago- 
nalization as a function of A for the case of 7 = 0.5. x 
shows different scaling behavior on two sides of A = 1. 
This is the quantum critical point of the LMG model. 
For A > 1, the system is in the polarized phase and spins 
are polarized along the z-direction. In this phase, x f° r 
different system sizes collapse to a single curve indicating 
that x ~ -^°- F° r A < 1, the model is in the symmetry- 
broken phase. The spin-spin interaction in the model in 
the xy plane dominates. The normalized x f° r differ- 
ent size collapse into a single curve and one has \ ~ N. 
These results are consistent with the scaling analysis us- 
ing the analytical expression of x as shown in Eq. (|38[) 
and Eq. P0"1) . respectively. 

To interpret the results obtained, we argue that in or- 
der to observe the QZE in the symmetry broken phase, 
the number of measurements Q has to be of comparable 
order of N. In a thermodynamic system, this condition is 
hardly realized. However, in the polarized phase where x 
is independent of the system size, the QZE is much more 
easier to be observed in comparing to that in the symme- 
try broken phase. Physically, this may be understood by 
the fact that in the case where A > 1, the ground state 
of the model is given by the configuration in which spins 
are already fully polarized along the external field. 

Moreover, from Fig. [13 we also see that x shows some 
abnormal behavior at the quantum critical point of the 
LMG model. Fig. [7]shows a plot of the first derivative of 
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the normalized \ with respect to A near the vicinity of the 
critical point. d{x/N)/d\ exhibits a minimum near A = 1 
and this minimum's amplitude increases as the system 
size increases. From the inset of the figure, one finds 
d(x/N)/d\ ~ jV°.324±o.ooi i n t he thermodynamic limit 
where N — > oo, d(x/N)/d\ is expected to be divergent 
at the critical point. 

Fig. [5] and [S] shows the numerical exact diagonalization 
result of x and d(x/N)/d\ as a function of A respectively 
for the case of 7 = 0. The scaling behavior of x obtained 
here is the same as the case for 7 = 0.5. From the inset 
of Fig. |jjj we also have d( x /N)/d\ ~ AT - 328 ± - 001 , which 
is consistent with the case for 7 = 0.5 up to two digits. 



SUMMARY 



the system sizes in order to observe the QZE, and sim- 
ilar conclusion is obtained for the case of the symmetry 
broken phase in the LMG model; however for the po- 
larized phase of the LMG model, the QZE can be easily 
observed via frequent measurements. Since we know that 
the size of the system in the thermodynamic limit is of 
order 10 23 , in reality it is almost impossible to have such 
a high frequency of projective measurement by the cur- 
rent technology. In this sense, we can safely argue that 
under certain circumstances, the QZE may breaks down 
in the thermodynamic limit. 

Furthermore, we also highlight the abnormal behavior 
of % as a function of the external driving parameter across 
the quantum critical point of the Ising model and the 
LMG model. The possibility of detecting QPT by QZE 
is raised and further exploration would be of interest. 



In this report, we have investigated the QZE from 
the viewpoint of condensed matter physics. We have 
obtained the scaling behavior of x m two analytically 
solvable systems, namely the one-dimensional transverse- 
field Ising model and the LMG model. We have found 
that in the Ising model, the frequency of the projective 
measurement should be of comparable order to that of 
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